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Intermittency in the Gledzer-Okhitani-Yamada (GOY) model of turbulence is explained in terms 
of collisions of coherent soliton-like structures with a random background issuing from the desin- 
tegration of their predecessors. This two-fluid picture is substantiated by the elucidation of local 
dynamical mechanisms leading to anomalous growth of coherent structures, their detection in true 
signals involving forcing and dissipation, and an investigation of their statistics. 


From a theoretical point of view, one of the most challenging features of fully developed turbulence is an interplay 
between random almost Gaussian background and coherent ordered structures responsible for deviations from Gaus¬ 
sian statistics. This duality explains in particular why field theoretical methods meet great difficulties in capturing 
intermittency effects directly from Navier-Stokes equations, even if considerable progress has been made recently in 
the related problem of random advection [I|,||. Although coherent structures were visualized as sheets or tubes of 
vorticity [|| in the case of 3D-incompressible turbulence, little is yet known about the way they form, their degree of 
stability and as a consequence their statistical relevance. 

Addressing such issues in the simpler context of the so-called shell models of turbulence may help one to figure 
out new mechanisms of intermittency, possibly at work in the complete Navier-Stokes dynamics. It was first realized 
by Siggia [Q that the one-dinrensional character of these models favours the formation of coherent soliton-like pulses, 
whose amplitude grows in a self-similar way, as they move from large to small scales. Many years later, Parisi Q, in 
an unpublished work, envisioned the turbulent medium forming in these systems as a gas of interacting “solitons”, 
with a continuous spectrum of scaling exponents leading to mutifractality. However, we have shown recently || that 
genuine self-similar solutions of the equations of motion in the inertial range display an unique scaling exponent (to 
be denoted below as zq), provided they are localized in k-space. When zo departs enough from the Komogorov value 
2/3, the scaling of large deviations observed in full simulations of the corresponding shell model (including energy 
injection at large scales and dissipation at small ones) is well accounted for by the properties of these ideal objects. 
But the resulting statistics, almost unifractal, is very different from realistic turbulence. 

In contrast to the previous situation, z o turns out to be quite close to 2/3 in the case of the Gledzer-Okhitani- 
Yamada (GOY) model 0, in the range of parameters where it reproduces very well the multiscaling properties of 
real turbulent flow ||. Taken literally, this result would suggest weak intermittency effects unless the hitherto ig¬ 
nored dressing of coherent structures by their interaction with the rest of the flow helps to produce more singular 
fluctuations. Acceleration of time scales downward the cascade, which is a salient feature of Navier-Stokes dynamics 
well captured by shell models, makes the collision between two coherent structures very unlikely. We are thus led 
rather naturally to a two-fluid picture, where coherent structures form in and propagate into a featureless random 
background. Somewhat paradoxically, the good transfer properties of the GOY model turn out to be a crucial ingre¬ 
dient in the story, allowing the persistence of a rich activity on all scales between the passage of two intense events. 
Randomness of the background comes primarily from the large separation of time scales between the integral and 
dissipative scales. In this Letter, we present various facts supporting this new physical picture of intermittency. 


Equations of motion for the GOY model in the inertial range may be cast in the form : 


db n 

dt 


Q 2 (1 - e)6;_ 2 6;_ 1 + - Q~ 2 b* n+1 b * n+2 , 


(i) 


where the complex variable b n = k n u n should be understood as the Fourier component of the gradient velocity field 
at wavenumber k n = Q n and the integer n runs from 0 to + 00 . Throughout this paper, usual values of parameters 
e = 0.5 and Q = 2 will be assumed. This choice yields values of scaling exponents C, p of the statistical moments (|it ra | p ) 
very close to those predicted by the She-Leveque formula [joj : 


C p = h 0 p + d 0 [l - 0 P ], 

with ho = 1/9, do = 2, and f3 = (2/3) 1 / 3 . Equations ®) admit formally self-similar solutions of the type : 


( 2 ) 
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( 3 ) 


b n (t) = Q nz g(v = Q nz (t-t*)). 


Provided the scaling function g vanishes for v —* — oo and goes to a constant for v —> 0, (Q) describes nothing but an 
invading front reaching the smallest scales in finite time and leaving in its trail a spectrum of slope z in logarithmic 
scale. In the ^-representation (the most relevant, dynamically), such a front appears like a soliton-like localized 
structure. For z = 2/3, Kolmogorov scaling is recovered. 

It is convenient to build from the b n ’s an infinite-dimensional vector b, whose squared norm (b, b) = \bn\ 2 

would be the enstrophy in real flow. The trend towards enstrophy blow-up can be cured by concentrating on the 
dynamics of the unit vector C = ^ and working with a desingularizing time variable t such that ^ = yj (b, b). 

By doing so, it was observed in |(]] that every initial condition of finite enstrophy leads to the same asymptotic state, 
up to a time translation, of the form b n (r) = e A ° T <5>o{n — vqt). Further, the envelope function <F 0 turns out to be 
purely real and positive up to trivial phase symmetries of the model |lo| leaving unaffected the energy flux e„ from 
shell to shell defined as : 


£n — Q 2n 3? [(1 — e)b n -2b n -ib n + b n -ib n b n +i\ . (4) 

From the easily measurable values of the velocity vo and the growth rate Aq, the estimate zo = q = 0.72 could 
be deduced. This is much lower than the exponent of the most intense event observed in full simulations of the GOY 
model, which, according to (|j), should be z max = 1 — ho = 0.89. 

To bring out the rules of interaction of such coherent pulses with a turbulent background, we first let one of them 
collide with some localized activity residing on a shell downstream. The system was prepared in a state consisting of 
a well-formed pulse (of unit norm and real positive) and, three shells in front of its center, a “defect” of amplitude 
a and relative phase 9 corresponding to an initial perturbation of the form Sb n = aexp(i0)(5„ ino . The dynamical 
rescaling method outlined before makes it very easy to compare the amplitudes of the pulse, with or without collision, 
when it crosses a shell far beyond the place of the collision. The relevant information is contained in the logarithmic 
amplitude gain G(a, 9) defined as : 

\b coll \ 

G(a, 9) = lim log . (5) 

B »"o K\ 

Fig. 0 shows the variation of G with 9 for two values of the amplitude, whose first one, a = 0.2, is still in the linear 
response regime. The first interesting observation is that a strong phase mismatch is required to enhance the growth of 
the pulse; the higher the strength of the collision, the more stringent gets this condition. Fig. ^ shows the behaviour 
of G as a increases, in the optimal case 9 = n. What ultimately limits the growth of the incident pulse is the fact 
that the new structure forming in front of it, as a result of the collision, gets too fast and starts to lead its own life. 
By contrast, we see in Fig. ||how for a = 1.5, i.e. slightly below the stability threshold, the latter is finally caught up 
by the former. 

In order to induce corrections in scaling exponent z, such collisions should occur repeatedly all along the cascade. 
It is reasonable to think that collisions with “defects” of large amplitude will keep their high efficiency only if they 
are sparse enough, because they involve long-lived intermediate states presumably prone to splitting instabilities. To 
get a better feeling for the order of magnitude of exponent corrections A z one may expect from such a mechanism, 
we just divided G(a, 7r) by the number A n(a) of steps the original pulse has to go forward before recovering its initial 
shape. This conservative estimate yields corrections in scaling exponents which can get as high as 0.16 or 0.18 on the 
An = 2 and An = 3 branches, corresponding respectively to 0.2 < a < 0.6 and 0.6 < a < 1.2. Those are the orders 
of magnitude expected to bridge the gap between zo and 1 — hq. 

The analysis of true signals corroborates to a large extent the two-fluid picture we propose. In order to detect 
coherent structures, we just selected series of local maxima of e n , starting from the top of the cascade and going 
downwards, such that at each step \b n \ grows by a factor larger than some prescribed value, conveniently written as 
Q z . By letting the effective exponent z vary, one scans events of various singular strength. Furthermore, coherency is 
controlled by demanding that the time interval t n +i — t n between the occurence of local extrema on neighbouring shells 
never exceeds the local turn-over time ~ Q~ nz compatible with the scaling of the singular event under examination. 
Figure 0 depicts such an event at three successive times for a Reynolds number Re = 10 6 . The imposed growth factor 
was Q °^ 5 , which this particular pulse failed to achieve well before reaching the dissipative shell of index rid = 15. The 
left side of Fig. showing \b n \ versus n, reveals that the coherent structure emerges from a disorganized K41-ramp, 
whose level does not seem to vary significantly during the time of observation (the length of the cascade is however 
too short to make a definitive statement on this last point, whose importance will be stressed later on). On the right 


2 







side of the same figure are shown the corresponding patterns of the gauge invariant phase ’F rl = 2 + 0 n -i + 0 n 

entering into the expression ( 0 ) of e n . One notes the reduction of the amplitude of spatial oscillations of the phase in 
the trail of the coherent structure, as well as the presence of a well-established phase defect just in front of it, able to 
induce a local change of sign in the energy flux. The same features were observed for every singular event we analyzed 
in this way. Altogether, they fit in nicely with the conclusions drawn from the study of elementary collisions. They 
also strongly suggest that the random background, necessary to feed anomalous growth of the coherent component, 
may simply obey mean-field Kolmogorov scaling. 

At this point we are facing an interesting dilemma. If incoherent fluctuations were obeying some intrinsic scaling 
fixed once for all, say of Kolmogorov-type, we would expect unifractal statistics in the limit of infinite Reynolds number 
because only the scaling exponent zq would survive to the increase of the cascade length. Multiscaling properties of 
the GOY model can be preserved asymptotically within our two-fluid picture, if and only if the turbulent background 
manages by a way or the other to stay at par with the coherent structure propagating into it. The simplest way to 
describe this situation is provided by the following stochastic dynamical system : 

= N[b] + rj , (6) 

where b now embodies the coherent part of the velocity field, N[b] is the nonlinear kernel of (0) and r/ is a Gaussian 
random force, delta-correlated in time, whose correlations read : 

(vUtfVn'it')) = T(b,b Y /2 5 nn '5(t-t'). (7) 

with T small for the sake of consistency. 

Within an adiabatic approximation keeping the shape of the soliton frozen, dynamics (0) reduces to the biased 
Brownian motion of two collective variables n(r) and B(t), respectively the position of the soliton and the logarithm 
of its amplitude (r is the proper time introduced before). The density of probability P n (z) for developing an effective 
growth exponent z after n steps appears then as a sum over all random walks starting from the origin (i.e. n(0) = 
B{ 0) = 0) and such that B(t) = nzlogQ at the time r of their first visit to the n th shell. For n> 1, P n {z) takes the 
form expected within multifractal descriptions |TT| : 

P n {z) ~ y/nexp[ns(z) + si(z)], (8) 


where the Cramer’s function s(z) reads : 


s(z) 


(z/z 0 ) +a 2 - \/((~/ Zo ) 2 + a 2 ) (1 + a 2 ) 


(9) 


In the above equation, a 2 = and c = AoZologQ are two numbers related only to the properties of the 

ideal self-similar solution. The bad thing with (0) is that it predicts a linear decrease of s(;z) at large 2 which we 
believe to be an artefact of the approximation used to solve (|0) and, in any case, is not confirmed by numerics. 

We tried indeed to extract s(z) from true signals, by investigating the statistics of local maxima of e n exceeding 1 (a 

value of the order of the average energy flux). Scaling exponents z were computed as 2 = | + 3 ^- 2 ) 

where t ^ denotes the time of occurence of the i th maximum on the n th shell, and j(i) is the greatest integer such 
that all shells of index between 2 and n present at least one local maximum within the time interval [t^)yt[ n ^]. For 
Re = 10 6 , the resulting pdf P n - 2 (z) was observed to vary with n in a way compatible with (JsJ) only in a very narrow 
window of shell indices 7 < n < 9. We plot in Fig. 0 s (n)(,) _ j 0 g a_| for n = 8 and 9, which should give 

a fair account of the asymptotic Cramer’s function s(z). It is worth noting that deviations from a parabolic shape 
are hardly perceptible for values of z as high as 0.87, i.e. not very far from the maximal exponent postulated in 
She-Leveque’s approach. However, there is no hint of an abrupt cut-off on z, in contradiction with the conclusions of 
0. It would be interesting to see whether fully nonlinear treatments of the stochastic dynamical system (Q), based 
on instanton techniques similar to those recently introduced in the context of turbulence [Q , yield a better agreement 
with Fig. 0 than the analytical result (jj|) does. Work in this direction is in progress. 


) 


To conclude, let us emphasize that our theory, though still calling for improvement, rests on robust physical 
ingredients which could well have their counterparts in Navier-Stokes dynamics. The renormalization of coherent 
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structures by their environment turns out to be so strong that the existence of ideal self-similar objects, leading to 
finite time blow-up of the enstrophy in the zero-viscosity limit, is certainly not a prerequisite for the whole approach 
to make sense. It remains to see how the random background should be modelled in the context of Navier-Stokes 
equations. 
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FIG. 1. The logarithmic gain G(a,9) versus 9 for a = 0.2 and a = 1.0. 
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FIG. 2. The optimal gain G(a, n) versus the amplitude of the collision 



FIG. 3. Evolution of the normalized field C„(t) = 



during a collision of amplitude a = 1 


.5, close to the pulse-spliting 


threshold. 
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FIG. 4. Three succesive snapshots of the propagation of a pulse for Re = 10 6 . Amplitudes f>„ are depicted on the left side, 
together with power laws Q nz for z = 2/3 (solid line) and 2 = 0.85 (dashed line), while phases are shown on the right. 




FIG. 5. Estimates for the Cramer’s function 5 ( 2 ), extracted from pdf’s P n (z). Statistics was run over 3 x 10 4 turn-over times 
of large scales. 
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